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Abstract 

Based on the phenomenon of curvature-induced doping in graphene we propose a class of Peltier 
cooling devices, produced by geometrical effects, without gating. We show how a graphene nanorib- 
bon laid on an array of curved nano cylinders can be used to create a targeted and tunable cooling 
device. Using two different approaches, the Nonequlibrium Green's Function (NEGF) method and 
experimental inputs, we predict that the cooling power of such a device can approach the order of 
kW/cm 2 , on par with the best known techniques using standard superlattice structures. The struc- 
ture proposed here helps pave the way toward designing graphene electronics which use geometry 
rather than gating to control devices. 
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As electronics become smaller, one problem facing the industry is that of targeted, on- 
demand cooling, conventionally achieved by semiconductor-based Peltier coolers^ - -, con- 
taining arrays of n-type and p-type pellets. However, one limitation of Peltier coolers made 
with conventional fabrication techniques is that carrier concentration for each pellet is set 
by doping level and cannot be adjusted after fabrication, which in turn results in that the 
cooler's cooling power cannot be dynamically changed. Therefore, new materials and novel 
structures are desired to fabricate a cooler with dynamically changeable cooling power. 

Graphene, a two-dimensional atomic layer, is a potential candidate for a high-performance 
thermoelectric cooler. As a result of its unique energy spectrum, the doping level can be 
controlled by shifting the Fermi level either above (n-type) or below (p-type) the Dirac point. 
Therefore, remarkable efforts have been focused on both electrical^ - — and chemical^ 1 ^ meth- 
ods for fabricating graphene-based p-n junctions, which are the main constituents of a Peltier 
cooler. However, due to the geometry of these p-n junctions, the direction of heat transfer is 
parallel to the fabrication substrate. For this reason, the geometry of coolers has to be such 
that the p- and n-type regions are perpendicular to the surface which needs to be cooled, 
like a stack of books on a shelf. This extra level of manufacturing makes the devices less 
likely to be of use, especially in nanoscale applications. 

Whereas the currently available fabrication techniques of the graphene-based p-n junction 
are confined to surfaces, the three dimensional configurations of this flexible 2D membrane 
may eventually be exploited to develop integrated circuits with components not achievable 
before. In this paper, we show how a single, continuous graphene nanoribbon (GNR) may 
be used to create a nanoscale Peltier array without need for lithography or gating. The key 
ingredients in any Peltier cooler are an array of p-n junctions which are electrically in series 
yet thermally in parallel, such that the junctions which evolve heat are on one side of the 
device, and the junctions which absorb heat are on the other. Upon application of a current 
from left to right through the GNR shown in Fig. [TJ heat will be pumped from the junctions 
near the bottom of the structure to the junctions near the top. From two different methods 
which agree, we estimate the cooling power of the proposed device to be on the order of 
kW/cm 2 . 

Curvature- induced doping provides a route to creating the p-n junctions required by 
Peltier coolers. When the local curvature of graphene resembles that of a sphere (large 
mean curvature), or that of a saddle (small mean curvature), the change in orbital energies 
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and overlap integrals associated with the curvature causes the energy spectrums and the 
Dirac points to shift differently in these two cases, leading to n-type or p-type doping^. By 
draping a single Armchair metallic GNR over a curved cylindrical protrusion (such as a bent 
nanotube may provide, see Fig. E ffib) ), a region of large mean curvature is created next to a 



region of small mean curvature, and a p-n junction is created. 

In Fig. [TJ blue regions denote areas of large mean curvature and red regions denote areas 
of small mean curvature. The curvature-induced Dirac-point shift is given by^ 

$(x) = -3aa 2 (#(x)) 2 . (1) 

H(x) is the mean curvature of the surface at the point x 



2\r l r-i, 

where t\ and r 2 are the principal radii of curvature at the point x, a is the nearest-neighbor 



distance, a = 2.5A, and a ~ 9.23eV. In our proposed device, see Fig. Ma) r\ = R% and 



r 2 = r for the outer surface while r\ = R\ and r 2 = r for the inner surface. Because the 
Dirac-point shift depends on the square of H~, it is independent of the coordinate system. 
The Dirac points in the regions with a larger magnitude of the mean curvature (locally like 
a sphere) are lowered more while the Dirac points in the regions with smaller magnitude 
of the mean curvature (regions which are flat, or alternatively which are saddle-like) are 
lowered less. This effect implies that regions with large magnitude of the mean curvature 
are locally electron-doped, while the regions with small magnitude of the mean curvature 
are locally hole-doped. 

The optimal cooling power P m for one Peltier cooling element can be derived from the 
rate of heat absorption 16 

P m = K^-ZTl - AT) = _ KAT, (3) 

where AT = T# — Tc is the temperature difference between the two heat reservoirs, AS* = 
S p — S n , is the difference between Seebeck coefficients in the p-doped (with smaller mean 
curvature) and n-doped (with larger mean curvature) regions, and Z = , the figure of 
merit of the p-n junction. K and R are the thermal conductance and electric resistance, 
respectively, of the p-n junction. When electrical current equal to / = I opt = ^ AS ^ Tc is driven 
through one Peltier cooling element, the optimal cooling power P m is reached. 
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FIG. 1. A single GNR laid on an array of curved nanotubes. A GNR bent into the geometry 
shown becomes a Peltier cooling device with cooling power «lkW/cm 2 . The large mean curvature 
of the blue regions leads to spontaneous n-type doping, and the low mean curvature of the red 
regions leads to spontaneous p-type doping. Upon application of a current from left to right, heat 
is pumped from bottom to top. The cooling power may be adjusted by changing the curvature, 
via, e.g., the application of uniaxial pressure on the device. 

In the following, we use two approaches to estimate the cooling power of our proposed 
device under the condition that AT = Tjj — Tc = 0. For the first approach, we use Nonequi- 
librium Green's Function (NEGF) method to perform calculations from the atomic scale. 
For the second approach, we estimate the cooling power from experimental measurements 
of the Seebeck coefficient^. Results from these two approaches are similar. 

The NEGF method can be used to calculate the transmission coefficient T (E) for a 
perfect GNR from the atomic scaled. The system contains one channel and two contacts. 
The Green's function G for the transport channel is denned as: 

G= [(E + i^I-H-Hx -X 2 ]-\ (4) 

where E is the incident energy of electrons, / is the identity matrix with the same 
dimension as that of the channel Hamiltonian H, rj is a small number accounting for the 
energy level broadening, and £i and S 2 represent the effects from two contacts. In order to 
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(a)two curvature radii (b)Procedures for building the proposed device. 

FIG. 2. (a). Creating two curvatures by bending. Two curvature radii with equal signs results 
in large mean curvature while two radii with different signs results in small mean curvature. i?i 
and r are the two radii for the inner side and R2 and r are for the outer side of the tube. (b). A 
schematic illustration of making our proposed nanocooler. At first, a GNR is put on the surface of 
a insulating nanotube. Secondly, curvatures are introduced on the GNR by bending the nanotube. 
The curving angle 6 X (9 y ) defined by L x and R\ (L y and r). 

describe the tunnelling processes between the channel and two contacts, two functions, I\ 
and T 2 , are introduced for two contacts respectively, 

r 1 = i[E 1 -E+] ; r 2 = i[s 2 -s+]. (5) 

The transmission coefficient for a perfect GNR is then given as^ 

T (E) = Trace[T 1 GT 2 G + ]. (6) 

In real systems, scattering processes have to be taken into account. The electron-electron 
interaction usually causes phase relaxation, which destroys the coherence of the electronic 
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wave functions. The electron-phonon interaction and the electron-impurity interaction result 
in backscattering and momentum relaxation processes. For high voltage bias between two 
contacts, effects from inelastic scatterings are also involved. Dominant scattering processes 
are determined by the physical system under consideration. 

Under usual experimental conditions, where the voltage bias is small and the whole system 
is in the linear response regime, there is a way to take elastic scattering effects into account 
within NEGF. For quantum transport, it is well known that the transmission coefficient 
T(E) can be related to its ballistic limit counterpart T (E) in the following way^, 

T(E) = j^ToiE), (7) 

where L is the system size along the transport direction and A is the backscattering mean-free 
path and is assumed to depend only on disorder and not on the incident energy E. Therefore, 
once Tq(E) is obtained, it is straightforward to calculate T(E). For this paper, calculations of 
the the transmission coefficient are based on the GNR system with L x = 75nm, L y = 25nm, 
and A = 400nm. 

After obtaining T(E) we can calculate the Seebeck coefficient S and the electrical con- 
ductance g. An intermediate function L n is defined as^ 

Lnfr, T) = U dET(E)(E - ^n- df{E d ^ T) ), (8) 

where h is the Plank constant, and / is the Fermi distribution function. Based on these 
intermediate functions, g and S can be computed as 

g(fi) =e 2 L (/i,T), (9) 

1 Li(/i,T) 

s ^ = 7tT^t Y (10) 

where e is the electron charge, \x is the chemical potential and T is the average temperature. 
In FigJ31 we show the electric conductance and the Seebeck coefficient at T = 300-fT as 
functions of the on-site voltage for clean and disorder cases. We find that the electrical 
conductance of the clean case is slightly larger than that of the disordered case while the 
Seebeck coefficient behaves similarly in both cases. We furthermore find that the maximum 



of the Seebeck coefficient is on the order of lOOfiV/K, consistent with the experimental 
measurements for Graphene^. 

The cooling power P of our proposed device can be evaluated using the calculated Seebeck 
coefficient and electrical conductance with disorder. When the surface of a GNR is curved, 
Dirac points for different regions are shifted up or down, similar to applying a gate voltage. 
Therefore, the calculated Seebeck coefficient as a function of on-site voltage can be used as 
the Seebeck coefficients for different Dirac-point shifts caused by the curvatures. Combined 
with the calculated electric conductance, the cooling power can be calculated by using Eqj3] 
with AT = 0. 

In Fig. HI we show the cooling power as a function of two curvature radii R = R\ and 
r, where i? = 4 £ ,r = -^ £ ,T = 3007^, and the distance between two nano tubes d = 15nm. 

1 Ox 1 vy 

L x and L y are fixed and < 9 x ,9 y < ir. Therefore R and r have lower limits ^* and — , 
respectively. As we expect, the maximum cooling power takes place for small R and r, 
which is on the order of kW/cm 2 . As an example, if we take as inputs L x = 75nm = 9 X R, 
L y = 25nm = nr, and 9 X = |, then the corresponding Dirac-point shifts for the inner side 
($i) and outer side ($2) can be obtained after considering charge neutrality : $! « 1.96 
mV, and $2 ~ —1-96 mV. By using the calculated Seebeck coefficient and the electrical 
conductance (Fig. [ ^(d) and Fig. [ ^b) respectively), we find that the cooling power is 



0.3kW/cm 2 . When the cooling device is curved more, 9 X = the cooling power can be 
estimated in a similar way to be 0.5kW/cm 2 . This shows that the cooling power can be 
tuned by changing the curvature of the device, for example, by applying uniaxial pressure. 

We now use an entirely different approach to estimate the cooling power of our proposed 
device, in which the Seebeck coefficients are taken directly from experimental measure- 
ments(Fig. 3(b) in ). Given the applied gate voltage V g , the resulting Dirac-point shift of 
the Graphene sample V r can be derived as^ 

where t is the thickness of the S1O2 substrate, vp is the Fermi velocity of graphene, and eo 
and e are the permittivities of free space and S1O2, respectively. By using experimentally 
determined Seebeck coefficient-^ and electrical conductivity—, the cooling power can be 
similarly estimated from Eq. |3j 

In Fig. |5l we report the cooling power as functions of the two curvature radii r and R as 
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(a)The electric conductance for perfect GNR. 
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(b)The electric conductance for disordered GNR. 
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(c)The Seebeck coefficient for perfect GNR. (d)The Seebeck coefficient for disordered GNR. 

FIG. 3. The sample size is that L x = 75nm, L y = 25nm and T = 300/f.(a). The electric 
conductance for the perfect GNR(b). The electric conductance for the disordered GNR with 
A = 400nm.(c). The Seebeck coefficient as a function of gate voltage for the perfect GNR.(d). The 
Seebeck coefficient as a function of gate voltage for the disordered GNR with A = 400nm. 



we do for the first approach. Results are consistent with that from the first approach, which 
is on the order of kW/cm 2 . As an example, if we take as inputs L x = 75nm = 6 X R, L y = 
25nm = irr, and 9 X = |, then the corresponding Dirac-point shifts for the inner side ($i) 
and outer side ($2) can be obtained: $1 ~ 1.96 mV, and $2 ~ —1-96 mV. This corresponds 
to applied gate voltages of $i 5 ~ 1.38 V, and $ 2s ~ —1.38 V. Given these two gate voltages, 
the Seebeck coefficients in both the n-type and p-type regions^ can be obtained. Combined 
with experimental data for the electric conductance ( lO 6 ^^ , corresponding to the mobility 
the cooling power P can be estimated to be P w When the cooling 
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device is curved more, $ x — from similar calculations, the cooling power is obtained to 




FIG. 4. The cooling power as a function of two curvature radii for the GNR with L x = 75nm, 
L y = 25nm, and A = 400nm, which is from NEGF calculations. 

be « 0.9-^. This result again shows that we can tune the cooling power by bending the 
nanotubes through, e.g., applying uniaxial pressure. 

In conclusion, we have proposed a graphene-based nano mechanical cooling device and 
estimated its cooling power using two different approaches: the NEGF method and exper- 
imental inputs. As a result of geometry alone, a series of P-N junctions are created in the 
proposed device such that by applying electric current, heat can be pumped perpendicular 
to the surface of the substrate. We find P ~ 0.5-^r, close to that achievable with the best 

cm* ' 

cooling devices ~ 1-^^. Most importantly, the cooling power of the proposed device can 
be adjusted by changing the curvatures, via, e.g., applying uniaxial pressure to the device. 
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FIG. 5. The cooling power as a function of two curvature radii, R and r for a GNR with L x = 75nm, 
L y = 25nm. The calculation is performed by using the experimental input^ 
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